Preprint typeset in JHEP style. - PAPER VERSION DCTP/01/30, IPPP/01/15, hep-ph/0104178 



One-loop QCD corrections to gluon-gluon 
scattering at NNLO* 



E. W. N. Glover" and M. E. Tejeda-Yeomans" 

a Department of Physics, University of Durham, Durham DH1 3LE, England 
E-mail: E.W.N. GloverOdurham . ac . uk, M . E . Te j eda-Yeomans@durham .ac.uk 



Abstract: We present the O (af) virtual QCD corrections to gluon-gluon scattering 
due to the self-interference of the one-loop amplitude. We give analytic expressions 
renormalised in the MS scheme and work in conventional dimensional regularisation. 
We write the structure of the infrared divergences from direct Feynman diagram 
evaluation in terms of the Catani formalism for infrared divergences. Formulae for 
the finite remainder are given in terms of logarithms that are real in the physical 
region. These results, together with those previously obtained for quark-quark and 
quark-gluon scattering complete the one-loop matrix elements for 2 — > 2 processes 
needed for the next-to-next-to-leading order contribution to inclusive jet production 
at hadron colliders. 

Keywords: QCD, Jets, LEP HERA and SLC Physics, NLO and NNLO 
Computations. 



*Work supported in part by the UK Particle Physics and Astronomy Research Council and by 
the EU Fourth Framework Programme 'Training and Mobility of Researchers', Network 'Quantum 
Chromodynamics and the Deep Structure of Elementary Particles', contract FMRX-CT98-0194 
(DG 12 - MIHT). M.E.T. acknowledges financial support from CONACyT and the CVCP. We 
thank the British Council and German Academic Exchange Service for support under ARC project 
1050. 



1. Introduction 



To improve our understanding of jet production in current and future high energy 
collider experiments we need more accurate perturbative calculations in quantum 
chromodynamics (QCD). The progress of recent years means that nowadays a Next- 
to-Leading Order (NLO) perturbative calculation is standard (see for example [1,2]). 
Developments of efficient techniques for the evaluation of matrix elements and nu- 
merical algorithms that allow the virtual and bremstrahlung contributions to be 
combined have made NLO perturbative calculations commonplace. However, new 
improvements in detector technology and an increase on statistics and luminosity re- 
quire more precise theoretical calculations. The reduction of theoretical uncertainties 
is crucial to refine our understanding of QCD and to identify signatures of possible 
new physics beyond the Standard Model. 

The calculation of fixed next-to-next-to-leading order (NNLO) virtual corrections 
for massless 2^2 processes is now an attainable task since a vital breakthough 
has been made in the calculation of two-loop master integrals [3, 4, 5, 6, 7, 8, 9, 
10, 11], in methods for reducing tensor integrals to master integrals [5, 6, 9] and 
in the prediction of the structure of infrared divergences [12]. The first two-loop 
contributions for light-like 2—^2 scattering were addressed by Bern, Dixon together 
with Kosower [13] and Ghinculov [14] who have calculated the maximal-helicity- 
violating two-loop amplitude for gg — > gg, and the QED processes e + e~ — > 
and e + e~~ — > e~e + , respectively. More recently, Gehrmann and Remiddi [15] have 
provided analytic expressions for the complete set of two-loop integrals with one off- 
shell leg. These integrals are sufficient to allow the two-loop amplitudes for three jet 
production in e + e~ annihilation to be evaluated. 

Recently, we provided analytical expressions for the O (a*) one-loop and two- 
loop contribution to quark-quark scattering [16, 17, 18], and quark-gluon scatter- 
ing [19], as well as crossed and time- reversed processes, in the limit where the quark 
mass is neglected. The O (a*) two-loop contribution to gluon scattering is given in 
Ref. [20]. To complete the list of ingredients needed for the NNLO virtual corrections 
to gluon-gluon scattering, the interference term of one-loop amplitude with one-loop 
amplitude needs to be included. In this paper we present analytical formulae for this 
contribution using conventional dimensional regularisation (CDR) and renormalised 
in the MS scheme. 

The organisation of our paper is the following. We provide our notation and 
describe our method briefly in Sec. 2. In Section 3 we give analytic expressions 
for the interference of one-loop amplitude with one-loop amplitude. In Sec. 3.1 
we follow closely the notation used in Refs. [12, 20] to write the infrared singularity 
structure in the MS scheme for the one-loop contribution to the NNLO calculation in 
terms of the one-loop bubble and the one-loop box integral given in Appendix A. We 
show that the anticipated singularity structure agrees with our explicit calculation. 



1 



Finally, in Sec. 3.2 we give the finite contribution in terms of logarithms that have 
no imaginary parts. We conclude with a brief summary of results. 



2. Notation 

As in Ref. [20] we consider the scattering process 

g(pi) + g(p2) + g(p 3 ) + gfa) -> o, (2.1) 

where the gluons are all incoming with light-like momenta and carry colour indexes, 
Oj, in the adjoint representation. 

Their total momentum is conserved and satisfies 

Pi +P2 +P3 + Pa = 0, Pi=0, 

where the associated Mandelstam variables are given by 

s = (pi+p 2 ) 2 , t=(p 2 +p 3 ) 2 , u =( Pl +p 3 ) 2 , s + t + u = 0. (2.2) 

We treat all external gluon states in D dimensions, i.e. we use CDR, and we 
renormalise the ultraviolet divergences in the MS scheme. In this scheme, the renor- 
malised four point amplitude can be written as 

\M) = 4™ s \\M<®) + \M^) + g) 2 \M^) + O (afj , 

(2.3) 

where a s = a s (fi 2 ) is the running coupling at renormalisation scale fi and the 
represents the colour-space vector describing the renormalised i-loop amplitude. The 
dependence on both renormalisation scale /j, and renormalisation scheme is implicit. 
The squared amplitude summed over spins and colours is denoted as 

V(s, t,u) = (M\M) = ^2\M(g + 9 -> 9 + g)\\ (2.4) 

which is symmetric under the exchange of s, t and u. 

Following from Eq. (2.3), the function V can be expanded perturbatively to yield 

V\s, t, u) + (f\ X? 6 ( S , t, u) + fc) 2 V*(s, t,u) + 



V(s,t,u) = 16n 2 a 2 s 



(2.5) 



where 



V\s,t,u) = (M^\M (0) ) 



16yiV 2 (l- e ) 2 (3-^-^-4V (2.6) 



s 



t 2 U 2 J 



V 6 {s,t,u) = (M {0) \M {1) ) + (.M (1) |.M (0) >, (2.7) 
V s (s,t,u) = (M^IM^) + (.M (0) |.M (2) > + (M {2) \M {0) ), (2.8) 
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and V = N 2 — 1 with N the number of colours. Expressions for T> % are given 
in Ref. [21] using CDR to isolate the infrared and ultraviolet singularities. Also, 
analytic formulae for the two-loop contribution to T> 8 

V^ 2x0 \s,t,u) = (M<®\MW) + (M^\M^) (2.9) 

are given in Ref. [20]. 

We now present expressions for the infrared singular and finite contributions to 
V 8 due to the interference of one-loop amplitudes with one-loop amplitudes, namely 

V 8 ^ lxl \s,t,u) = (M (l) \M (l) ). (2.10) 

This contribution is somewhat simpler to evaluate than the two-loop contribution 
but it is a vital part of the NNLO virtual corrections and we present it here for 
completeness. 

As in Refs. [16, 17, 18, 19, 20], we use QGRAF [22] to produce the one-loop Feynman 
diagrams to construct I.M*- 1 -'). We perform the internal trace (for fermion loops) 
over Dirac matrices in D dimensions. The tensor integrals that emerge from this 
amplitude are associated with scalar integrals in higher dimensions and with higher 
powers of propagators [9, 23]. 

These integrals, in turn, can be reduced to master integrals in D = 4 — 2e with 
the systematic application of integration-by-parts identities [24]. At the end, all 
integrals arising in the one-loop amplitudes can be written in terms of the one-loop 
bubble integral in D = 4 — 2e and the finite one-loop box integral in D = 6 — 2e 
(see Appendix for expansions around e = 0). This choice of master integrals is not 
unique but is useful to have a natural separation of the infrared poles and the finite 
part of the one-loop amplitudes. 

Finally, we project by (.M^l and perform the summation over colours and spins. 
Note that when we sum over the gluon polarisations, we ensure that the polarisation 
states are transversal by the use of an axial gauge 

]T efer = -g^ + ^ + U ' P ' (2.11) 

spins U i ' Pi 

where pi is the momentum, e« is the polarisation vector and n« is an arbitrary light- 
like 4-vector for gluon i. For simplicity, we choose = p%, n 1 ^ = pi, n% = p^ and 

n 4 — F3 ■ 

3. One-loop contribution 

We further decompose the one-loop contributions as a sum of two terms 

V 8(lxl) (s,t,u) = Voles(s,t,u) + Finite(s,t,u). (3.1) 
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Voles contains infrared singularities that will be analytically canceled by those oc- 
curring in radiative processes of the same order (ultraviolet divergences are removed 
by renormalisation). Finite is the remainder which is finite as e — > 0. 

3.1 Infrared Pole Structure 

Following the procedure outlined in Refs. [12, 20], we can write the infrared pole 
structure of the one-loop contributions renormalised in the MS scheme as 



Voles = 2 Re 



_I(^(o)| i -d)t (e)J (i) (e) |_ M (o) ) 



(M^\I^(e)\M^) + (M^\I^(e)\M^) 



(3.2) 



where |A^^) is the tree amplitude and \M^ 1,un ^) is the unrenormalised one-loop 
amplitude, obtained by direct Feynman diagram evaluation. Also, the first coefficient 
of the QCD beta function for Np (massless) quark flavours is 



11C A -AT R N F 
6 



with 



T, 



R 



1 

2" 



(3.3) 



It is convenient to decompose \M^) and \M.( 1 > un )) in terms of SU(N) matrices in 
the fundamental representation, T a , [25, 26, 27, 28, 29] so that \M {0) ) and \M {1 > un) ) 
may be expressed as nine-dimensional vectors in colour space 



\M {1,un) ) 



(Ti, r 2 , r 3 , r 4 , r 5 , r 6 , o, o, o) T , 

(Ci, £2, £3, £4, £5, £&, £7, £8, A 



(3.4) 
(3.5) 



where T indicates the transpose vector. Here the tree and loop amplitudes % and 
£i are the components of \M.^) and lA^ 1 '*"^) in the colour space spanned by the 
(non-orthogonal) basis 

— Tr (r^T^T^T"" 1 ) 
_ ^ ^2~ ,ai X' a2 T a4 T a3 ) 
_ ^ ^2~ ,ai X' a4 T a2 T a3 ) 

_ fj^j, ^2~" ll 2~ ,a3 2~ ,a4 2~ ,a2 ) 

— r J^j. /^pa\rpa4rpairpa2^ 

C 7 = Tr (T ai T a2 ) Tr (T a3 T a4 ) , 
C 8 = Tr (T ai T a3 ) Tr (T a2 T a4 ) , 

C 9 = Tr (T ai T a4 ) Tr (T a2 T a:i ) . (3.6) 
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Note that the amplitudes themselves are not required since we compute the interfer- 
ence of tree and loop amplitudes directly. 

In the same colour basis, the infrared-singularity operator J^(e) introduced by 
Catani [12] has the form [20] 



(T - U) (S - U) \ 
(U - T) (S - T) 
(T — S) (U — S) 
(T - S) (U - S) 
(S - T) 
(S - U) 


2MJ 








2N1 
(3.7) 



where 



S = 



s + iO 




T= 




U= 



(3.8) 



The matrix 1^ (e) acts directly as a rotation matrix on \M {0) ) and \M {1 ' un) ) in colour 
space. Note that on expanding S, imaginary parts are generated, the sign of which 
is fixed by the small imaginary part -MO. Other combinations such as (e) 
are obtained using the hermitian conjugate operator J^^(e) where the only practical 
change is that the sign of the imaginary part of S is reversed. 

The contraction of the colour vector \X) with the conjugate colour vector (Y\ 
obeys the rule 

(Y\X) = E E E YrXjC*C r (3.9) 

spins colours i,j=l 

In evaluating these contractions, we typically encounter I] co i ours C* Cj which is given 
by the ij component of the symmetric matrix CC 

^ C\ Ci Ci Ci C*2 

c 3 
c 2 
c 2 

Ci 

c 2 



00 



V 
16N 2 



C 2 
C 2 
C 2 
C 2 

c 3 



c 2 

Ci 

c 2 
c 2 
c 3 
c 2 



c 2 
c 2 

Ci 

c 3 
c 2 
c 2 



c 2 
c 2 
c 3 

Ci 

c 2 
c 2 



c 3 
c 2 
c 2 
c 2 
c 2 



-N 
NV 



NV 
NV 
-N 
-N 
NV NV 



NV \ 
-N 



NV NV 
NV NV 



NV 



NV NV -N -N NV NV N 2 V 
-N NV NV NV NV -N N 2 
NV -N NV NV -N NV N 2 



-N 

N 2 
N 2 V 
N 2 



-N 
NV 
N 2 
N 2 
N 2 V 



(3.10) 
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with 

C\ = iV 4 — 3iV 2 + 3, C 2 = 3-N 2 , C 3 = 3 + N 2 . (3.11) 

Similarly, we find that the interference of the tree-level amplitudes Sspins %*% ls 
given by TTy, where 



TT = 6^l-e) 2 {t 2 + ut + u 2 ) 2 



(3.12) 



s 2 t 2 u 2 
and the vector V is 

V = (u, t, s, s, t, u, 0, 0, 0). (3.13) 

Also, the interference of the tree-level amplitudes with one-loop amplitudes Lupins 
is given by TLij, where 

TC = V T W, (3.14) 

and the vector W is 

W = (F(s,t), F{s,u), F{u,t), F{u,t), F(s,u), F(s,t), G, G, G). (3.15) 

The function JF(s, t) is symmetric under the exchange of s and t, while G is symmetric 
under the exchange of any two Mandelstam invariants, so that 

F(s,t) = f 1 (s,t,u)+f 1 (t,s,u), (3.16) 
G = f2(s, t, u) + f 2 (s, u, t) + f 2 (t, s, u) + f 2 (t, u, s) + f 2 (u, s, t) + f 2 (u, t, s). 

(3.17) 

Here f\ and f 2 are given in terms of the one-loop box integral in D = 6— 2e dimensions 
and the one-loop bubble graph in D = 4 — 2e (see Appendix A), 

f( , , 167V(l-2e) 



+ 



+ 



f2(s,t,u) 



L 2(l - e) 2 (s 4 + s 3 t + st 3 + t 4 ) + 3(1 - 5e)s 2 t 2 ] Box 6 (s,t) 

8Nf{1 - 2e) [(1 - e) 2 (s 2 + t 2 ) + e(l + 3e)st] Box 6 (s, i) 
st 

16iV(l - 6) r , l2 _ + + s 4 + / 24 _ 5ge + 5Qe 2 _ ^3 _ g3f 

s 2 t 2 ue(3 — 2e) LV 7 v 7 

+ (36 - 99e + 93e 2 - 24e 3 - 2e 4 ) s 2 t 2 + (1 - e) (24 - 50e + 23e 2 ) st 3 

+4(1 - e)(l - 2e)(3 - 2e)t 4 ] Bub(t) 
16iV F 
st 2 u(3 - 2e) 

+ (6 - 15e + 21e 2 - 8e 3 ) st 2 + (1 - e) (5 - 6e + 2e 2 ) t 3 ] Bub(t), (3.18) 
-4(1 - efst + 3(1 - 5e)u 2 } Box 6 (m, i) 



(4 - 12e + 16e 2 - 4e 3 ) s 3 + (3 - lOe + 23e 2 - 8e 3 ) s 2 t 



32(1 -2e) 



32(1 -e) 



esu A 



4(1 - 2e)(l - e)t 2 + (8 - 17e)(l - e)ut 



+ (6 - 20e + 15e 2 + e 3 ) u 2 ] Bub(s). 



(3.19) 



6 



It can be easily noted that the leading infrared singularity in Eq. (3.2) is O (1/e 4 ). 
It is a very stringent check on the reliability of our calculation that the pole structure 
obtained by computing the Feynman diagrams directly and introducing series expan- 
sions in e for the scalar master integrals agrees with Eq. (3.2) through to O (1/e). 
We therefore construct the finite remainder by subtracting Eq. (3.2) from the full 
result. 

3.2 Finite contributions 

The finite two-loop contribution to V 8 (s,t,u) is defined as 

Finite(s,t,u) = V 8(lxl \s,t,u) - Voles(s,t,u), (3.20) 

where we subtract the series expansions of both £> 8(lxl )(s, t, u) and Voles(s ) t ) u) and 
set e — > 0. 

For convenience, we introduce the following logarithms 

X = log (—), Y = log (—) , S = log ( -L), (3.21) 



\ s / V s J \/i 2 j 

where /x is the renormalisation scale. 

We choose to present our results by grouping terms according to the power of 
the number of colours iV and the number of light quarks N F , so that 

/ N 2 \ 

Finite(s, t,u) = Vl N A A + N 2 B + N 3 N F C + NN F D + N 2 N F E + N F F + j^G J , 

(3.22) 

where 

A = | i (x 2 - 2 X Y + Y 2 + tt 2 ^ (x 2 - 2 X Y - 2 X + Y 2 + 2 Y + tt 2 ^ 

+ Lx A -AX*Y -— X 3 + 6X 2 Y 2 + 20X 2 Y-— X 2 S + 10 X 2 n 2 + — X 2 
V 3 3 9 

154 785 

-AXY 3 -20XY 2 -AXYn 2 -6XY + XS-lQXn 2 + X + Y A 

9 27 

+ 4 Y 3 --Y 2 S + 2Y 2 n 2 --Y 2 + —Y S +1QY n 2 + —Y + — S 2 
3 9 3 9 9 

2948 „ 4 2 9014\ t 2 
+ U X 1 + 12 A' 3 + i A' 2 F - " X 2 S + 16 A' 2 tt 2 - y x 2 + ^ A" y 2 

o y o o y o 

2 536^ 116 ^ 2 4 8 4 02 44 2 5896 _ 58 2 18028 \ t 

+ 16 F 2 tt 2 F H Y 7r -\ S 2 S tt 2 H S H tt 2 H 

9 3 9 3 27 9 81 J u 
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+ X 2 (x 2 + 47r 2 j^ 

+ 2X (x 3 + 2X 2 + 4X7t 2 + 4t/) ^ 

+ ( 7 X^-8X 3 Y -—X 3 + 12X 2 Y 2 + —X 2 Y-—X 2 S + 2AX 2 7T 2 + —X 2 

V 3 3 3 9 

-8xy 3 --xy 2 + -xy5-8iy, 2 + -xy + -i5--x, 2 

3 3 9 9 3 

2948 , 88 ,~ 2 536 ta2 88 ^ 2 484 o2 44 „ 2 

+ X + AY Y 2 S + 12 F 2 7T 2 y H Y n + S 2 S 7T 2 

27 3 9 3 9 3 

5896 „ „ 4 10 2 18028 \ t 2 

+ fl!x 4 -iox 3 r--x 3 + — x 2 r 2 + nx 2 r-— x 2 1 s + 29X 2 tt 2 

V 2 3 2 3 

-i^ + 22.YKS-5.YYV + ^A-r + Hl. YS + M A ^ 
3 9 9 3 

5309 v 605 o2 2 7667 „ 5 4 113 2 24533 \ 1 f 1 



(3.23) 



5 = 1 6 (x 2 - 2 X y + Y 2 + vr 2 ^) (x 2 - 2 X Y - 2 X + Y 2 + 2 F + tt 2 ^ 

+ ^72 X 4 - 120 X 3 y - 356 X 3 + 48 X 2 Y 2 + 580 X 2 Y + 156 X 2 tt 2 + X 2 

+24Xy 3 -404Xy 2 - 144Xyvr 2 - i^Xy- 392Xvr 2 - 112X-24y 4 

3 

\ t 2 

+180 y 3 - 12 y 2 tt 2 - 32 y 2 + 392 Y tt 2 + 112 Y + 12 vr 4 + 12 tt 2 — 

y s 2 

+ (^-24X 4 + 144X 3 y + 408X 3 - 48X 2 Y 2 - 272 X 2 Y + 120 X 2 tt 2 -64X 2 

2752 

+96 X y 3 + 624 X y 2 + 288 X y tt 2 + X y + 792 X tt 2 + 224 X 

o 

+144 y 2 tt 2 + 528 y tt 2 + — ^ L 

3 y u 

+12 X 2 (x 2 + 4vr 2 ^ 

+24X (x 3 + 2X 2 + 4Xvr 2 + 47r 2 j ^ 

+ ^84 X 4 - 96 X 3 y - 104 X 3 + 96 X 2 Y 2 + 352 X 2 Y + 288 X 2 tt 2 + X 2 
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1112 \ t 2 

-96 X Y tt 2 - 32 X tt 2 + 96 Y 2 n 2 + 352 Y n 2 + 24 tt 4 + tt 2 — 

3 y -u 2 

+ ^42 X 4 + 32 X 3 + 66 X 2 Y 2 + 41Q X 2 Y + 288 X 2 tt 2 + X 2 + 84 X Y tt 2 
+ — Xr + 716X7r 2 + 112X + 157r 4 + 6667r 2 + 48) l + ju^tl (3.24) 



| - |x 2 - 2 X Y + Y 2 + tt 2 ^X 2 - 2 X y - 2 X + Y 2 + 2 Y + tt 2 ^ 

+ f-x 4 + 2X 3 y + — x 3 -3X 2 y 2 -22X 2 y + — x 2 5-3X 2 tt 2 

\ 6 6 

185 x 2 + 2 x y 3 + 22 x y 2 + 2 x y tt 2 + 12 x y - — + - x tt 2 



18 "" ~ ~ ' 9 3 



27 2 2 18 3 3 9 

88 2 976 1 4 , 2752 \ t 2 

+ ^-4X 4 -4X 3 y-27X 3 + 6X 2 y 2 + 12X 2 y-X 2 S- 18X 2 tt 2 
127 9 „ , 9 148 62 

x 2 -4xy 3 -iixy 2 + 2xy 1 s-4xy7r 2 xy — xs 



3 9 3 9 9 

27 9 81 y u 

-2X 2 (x* + 4tt 2 )£ 

-4X ^X 3 + 2X 2 + 4Xvr 2 + 47r 2 ^ ^ 

+ ^-7X 4 -yX 3 -yX 2 y + ^X 2 1 S- 28X 2 7r 2 -yX 2 + ^Xy 2 

_^xy,-fxy-f x.-fx. 2 -f x + ^y 2 , + fy 2 

_16 yff2 _176 52+ 8 1952 40 5504^ 

3 9 3 27 3 81 J u 

+ ^-4X 4 + 2X 3 y-yX 3 -^x 2 y 2 -2X 2 y + ^x 2 1 s-i5X 2 7r 2 
131 x 2 -4xys + xy. 2 -^xy-^xs-^x. 2 -±^x 



18 9 9 3 27 

™ s i + 2S% >™ s _^_^- _ 9337)1 + L„ t l ( 3. 25) 



9 27 4 18 81 
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D = I - 2 |x 2 - 2 X Y + F 2 + TT 2 J |x 2 - 2 X F - 2 X + F 2 + 2 F + tt 2 J t — 

+ ( - 4 X 4 + 8 X 3 Y + 75 X 3 - 6 X 2 Y 2 - — X 2 F - 8 X 2 tt 2 - — X 2 

293 9 9 424 232 , 116 A , 

+— XF 2 + 4XFvr 2 + — XF + — Xtt 2 + — X + 2F 4 -43F 3 

o o o o 

+4F 2 7T 2 + — F 2 - — F 7T 2 - — F - 7T 4 - 4 TT 2 ^ ^ 

3 3 3 J s 2 

+ ( -4X 4 - 16X 3 F- 102X 3 + 12X 2 F 2 + — X 2 F-32X 2 vr 2 + — X 2 
V 3 3 

_ 8 XF 3 -— XF 2 -24XFtt 2 -312XF- — Xtt 2 -— X-4F 2 tt 2 
3 3 3 

332 „ 2 n 4 820 ,\ i 
Ftt 2 + 2tt 4 tt 2 ' 



3 3 J u 

-4X 2 (X 2 + 47 T 2 )^ 

-8X ^X 3 + 2X 2 + 4Xvr 2 + 47r 2 ^ ^ 

+ 14X 4 + 8X 3 -64X 2 F-56X 2 tt 2 - -y^ 2 - 16 X tt 2 - 64 F vr 2 

_4oo \ t 

3 y u 

( 35 
+ ( - 4 X 4 - 4 X 3 F - y X 3 - 85 X 2 F - 26 X 2 tt 2 - 206 X 2 - 2 X F tt 2 

1<S ir4i»' -¥x-U <-^-u>)UL»t) (3-26) 



3 

£?= |^ (x 2 -2XF + F 2 + 7r 2 ^J (x 2 -2XF-2X + F 2 + 2F + tt 2 ^ t — 

+ f-ix 4 + 2X 3 F-X 3 -3X 2 F 2 + 2X 2 F- I X 2 S - X 2 tt 2 + — X 2 

V 2 3 9 

10 2 22 1 

+2XF 3 -2XF 2 + 2XFtt 2 -6XF + — X5--X7T 2 X--F 4 

9 3 27 2 

+Y z + -Y 2 S-Y 2 tt 2 + -F 2 + -Y S+-Y n 2 + —Y +- S 2 + —S 
3 3 3 3 9 9 27 

1 4 2 236\ t 2 

7T + 7T H — 

2 si ; s 2 

+ fx 4 + 6X 3 --X 2 F + -X 2 5 + 4X 2 vr 2 + -X 2 + -XF 2 --XF 1 S 

V 3 3 3 3 3 

4. 4 vq 4. 34 ,68 4 2 16 2 2 , 160 
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22 2 472\ t 
3 81 J u 

+X 2 (^ + 4^^ 

+2X |x 3 + 2X 2 + 4Xvr 2 + 47r 2 ^ 

^+ ( 2 X 4 + 6X 3 + 8X 2 7r 2 + f X 2 + ^XS + 12X7r 2 + ^X + ^S 2 
w y y y ^ ( y 

160 n 44 2 472 \ t 2 

H S H vr 2 H — 

27 9 81 J u 

+ (}LX 4 -X 3 Y + -X 3 + -X 2 Y 2 -X 2 S + -X 2 n 2 + —X 2 --XYn 2 
\2 3 4 2 9 2 

9 9 3 27 9 27 8 

37 2 1049 \ If] 

+ T" + -ir r u "« (3 - 27) 



j - |x 2 - 2 X F + F 2 + tt 2 ^) |x 2 - 2 X F - 2 X + F 2 + 2 F + vr 2 ^ 

+ |x 4 - 4 X 3 F - 2 X 3 + 6 X 2 F 2 + y X 2 F + 2 X 2 vr 2 + 6 X 2 - 4 X F 3 

14 4 4 4 22 

XF 2 -4XF7r 2 + -XF--X7T 2 + -X + F 4 + 2F 3 + 2F 2 7r 2 F 2 

o o o o o 

+ f-2X 4 -4X 3 --X 2 F-8X 2 7r 2 --X 2 + -XF 2 + -XF--X7r 2 
V 3 3 3 3 3 

_!x + -F7r 2 + 167r 2 ^ - 
3 3 J u 

-4X ^X 3 + 2X 2 + 4Xvr 2 + 47r 2 ^ ^ 

+ ( - 4X 4 - 12X 3 - 16X 2 vr 2 - -X 2 - 24Xtt 2 + - A- 
\ 3 3 J u 2 

+ ( - X 4 + 2 X 3 F - - X 3 - - X 2 F 2 + 2 X 2 F - 3 X 2 vr 2 + — X 2 

+XF7T 2 + 2XF + -X7r 2 - -X- i7r 4 + — 7T 2 -8^) } + («<-> A (3.28) 
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G = 1 3 |x 2 - 2 X Y + Y 2 + tt 2 ^ (x 2 - 2 X Y - 2 X + Y 2 + 2 F + ^ ^ 

-3(X 2 -2XF + 2X + r 2 -2r + 7r 2 -2| (X 2 -2Xr-2X + F 2 + 2F + 7r 2 > t 



/ \ / - 2 

+ ^6X 4 + 24X 3 + 24X 2 7r 2 + 36X 2 + 48X7r 2 + 24X + 24vr 2 ^ ^ 

+ 12X ^X 3 + 2X 2 + 4Xvr 2 + 47r 2 ^ ^ 

+ ^12 X 4 + 36 X 3 + 48 X 2 7T 2 + 36 X 2 + 72 X 7T 2 + 24 tt 2 ^ ^ 

+ ^3 X 4 - 6 X 3 Y + 6 X 3 + ^ X 2 y 2 + 9 X 2 tt 2 + 6 X 2 - 3 X Y vr 2 + 6 X F 

+ 12Xtt 2 + 12X + 1 7T 4 + 6 7T 2 + 24^ | + |u <-> i| (3.29) 

We see that although we expect the finite piece to contain polylogarithms, they are all 
predicted by the infrared singular structure and are obtained by expanding Eq. (3.2) 
through to O (e). This is because the polylogarithms appear as the O (e) and O (e 2 ) 
terms in the expansion of the box integral in D = 6 and must be multiplied by an 
infrared singular term to contribute at O (1). At O (1), the interference of one box 
graph with another only collects the 0(1) terms in each and therefore yields only 
logarithms. 



4. Summary 

In this paper we presented analytic expressions for the O (a 4 ) QCD corrections to 
the 2^2 gluon-gluon scattering process due to the self-interference of the one-loop 
amplitude in the MS scheme. Throughout we employed conventional dimensional 
regularisation. 

The renormalised matrix elements are infrared divergent and contain poles down 
to O (1/e 4 ). The pole structure of the one-loop contribution is described by Eq. (3.2) 
using the formalism of Catani [12] while analytic formulae for the finite part according 
to the colour decomposition of Eq. (3.22) are given in Eqs. (3.23) to (3.29). The 
interference of the tree-level diagrams with the one-loop graphs are expressed in 
terms of the one-loop bubble graph in D = 4 — 2e dimensions and the one-loop box 
graph in D = 6 — 2e dimensions for which series expansions around e = are provided 
in Appendix A. 

The results presented here, together with those previously computed for quark- 
quark scattering [16, 18, 17], quark-gluon scattering [19] and gluon-gluon scatter- 
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ing [20] form a complete set of virtual matrix elements for parton-parton scattering 
at O (oig). They are vital ingredients for the next-to- next-to-leading order predictions 
for jet cross sections in hadron-hadron collisions. 

The next step is to combine these matrix elements with the tree-level 2 — > 
4 [30, 31, 32, 33], the one-loop 2^3 [34, 35, 36] contributions making sure that 
the infrared singularities are cancelled. This is a challenging task and a systematic 
procedure needs to be established. However, recent progress in determining the 
singular limits of tree-level matrix elements when two particles are unresolved [37, 38] 
and the soft and collinear limits of one-loop amplitudes [39, 29, 40], together with 
the analytic cancellation of the infrared singularities in the somewhat simpler case 
of e + e~ photon + jet at next-to-leading order [41], suggest that the technical 
problems may soon be solved for generic 2 — > 2 scattering processes. 

Initial state radiation complicates the issue and the colinear singularities from the 
incoming partons must be factorised into the parton density functions. The evolution 
of the parton density functions must also be known to an accuracy matching the 
hard scattering matrix element and requires knowledge of the three-loop splitting 
functions. Using the existing three-loop order moments [42, 43, 44], Van Neerven 
and Vogt have provided accurate parameterisations of the splitting functions in x- 
space [45, 46] which are now starting to be implemented in the global analyses [47]. 

Finally, and most importantly for phenomenological applications, an extension of 
the numerical next-to-leading order two jet programs [1, 2] to next-to-next-to-leading 
order must be developed. We note that Refs. [48, 49] have provided next-to-leading 
order programs for three jet production that form a natural starting point. We 
are therefore optimistic that these issues will shortly be resolved thereby enabling 
next-to-next-to-leading order QCD estimates of jet production in hadron collisions. 
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A. One-loop master integrals 

In this appendix, we list the expansions for the one-loop box integrals in D = 6 — 2e. 
We remain in the physical region s > 0, u, t < 0, and write coefficients in terms 
of logarithms and polylogarithms that are real in this domain. More precisely, we 
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use the notation of Eq. (3.21) to define the arguments of the logarithms and the 
polylogarithms are defined as 

r w dt. 



Li n (w) = / — Li n _i(t) for n = 2, 3, 4 
Jo t 



Li 2 (w) 



dt 

o 7 



log(l-t). 



(A.l) 



Using the standard polylogarithm identities [50], we retain the polylogarithms with 
arguments x, 1 — x and (x — l)/x, where 



x 



t u u x — 1 
y = = l-x, z = -- = . 

s s t x 

We find that the box integrals have the expansion 

\2 



(A.2) 



Box b (u,t) 



e^T(l + e)r(l-e) 2 ( tfY f 1 
2sr (1 - 2e) (1 - 2e) \ s J [2 

+2e 



Li 3 (x) - XLi 2 (x) - ^X 3 - *-X 



[x - yy + vr 2 

2 



-2e 2 



1 



\AAx) + FLi 3 (:r) - ix 2 Li 2 (x) - -X 4 - -X 3 Y + -X 2 Y 2 



^x 2 ff2 xv ^ 



+ («^t)j+0(e 3 ), (A.3) 



and 

Box 6 (s,t) 



e^r(l + 6)r(l-e) 2 
2ur(l-2e)(l-2e) \ u 



X 2 + 2mX) 



+e 



-2Li 3 (x) + 2XLi 2 (a;) - ^X 3 + 2YX 2 - tt 2 X + 2( 3 ) 



7T 



+irr 2Li 2 (x) + 4YX - X 2 



2U 4 (z) + 2Li 4 (y) - 2FLi 3 (a;) - 2XLi 3 (y) + (2XY - X 2 - 7r 2 )Li 2 (a;) 

+ Ix 4 - -X 3 Y + -X 2 F 2 + -tt 2 X 2 - 27r 2 XF + 2Y( 3 + -n 4 ) 
3 3 2 3 S3 6 J 

+m^-2Li 3 (a;) - 2Li 3 (y) + 2FLi 2 (:r) + ^X 3 - 2X 2 Y + 3XF 2 



7T" 



- Y + 20 



+ O (e 3 ) . 



(A.4) 



Box 6 (s,w) is obtained from Eq. (A.4) by exchanging u and t. 

Finally, the one-loop bubble integral in D = 4 — 2e dimensions is given by 



Bub(s) 



e^r(l + e)r(l- e ) 2 

r(2-2e)e \ s + i0 , 



(A.5) 
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